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(^ We propose a generalization of the two-dimensional eikonal-limit cloak derived from 

a conformal transformation to three dimensions. The proposed cloak is a spherical 



shell composed of only isotropic media; it operates in the transmission mode and 



\Q requires no mirror or ground plane. Unlike the well-known omnidirectional spherical 

cloaks, it may reduce visibility of an arbitrary object only for a very limited range 

Q of observation angles. In the short-wavelength limit, this cloaking structure restores 

"q^ not only the trajectories of incident rays, but also their phase, which is a necessary 

ingredient to complete invisibility. Both scalar-wave (acoustic) and transverse vector- 
O 



wave (electromagnetic) versions are presented. 



I. INTRODUCTION 

Transformation optics^ (TO) and transformation acoustics^'Sl (TA) are the two mathe- 
matically similar design methodologies based on the form-invariance of the vector and scalar 
Helmholtz equation, respectively. They have recently attracted significant attention as the 
tools that enabled novel ground-breaking applications, including manipulation of the appar- 
ent (effective) width of electromagnetic^^ (EM) and acoustic^"^ scatterers. The similarity 
between the equations describing EM and scalar acoustic waves is so close that many recipes 
derived for EM applications can be applied in the acoustics domain, and vice versa. In this 
article, we refer to transformation optics and acoustics collectively as transformation wave 
dynamics (TWD). 

Transformation wave dynamics, in its full generality, employs extremely general coordi- 
nate transformations, including such exotic options as non-bijective (space-folding)P^ and 
even complex-valued coordinate transformations'^^. However, the more general the choices 
of transformations are, the more complex the material properties required for their im- 
plementation tend to be. For example, non-bijective space-folding transformations require 
negative-index mediaf^, and complex-valued transformations typically lead to polarization- 
dependent gairi^. 

In light of this, it is reasonable to ask whether practically useful applications, such as EM 
or acoustic invisibility, can be accomplished with a more restrictive set of transformations, 
which lead to a narrower range of material properties. For example, one may wonder whether 
interesting TWD applications may arise from transformations that require only isotropic 
materials. The latter transformations are known to exist at least in two dimensions, where 
they form a subset of measure zero in the set of all possible real-valued transformations 
of TZ^. Indeed, conformal (angle-preserving) transformations are known to introduce no 
anisotropy into the material property tensors, a feature that has recently been explored for 
potential applications in radio engineeringHHSl, optics^ and acoustics^^"^. 

From the ease of fabrication perspective, one must notice that as long as a TWD de- 
vice requires exotic medium properties other than anisotropy, its implementation invariably 
relies on metamaterials^, photonic/phononic crystals^ or, generally speaking, artificial me- 
dia with naturally unavailable properties. In the case of optical cloaking, the other exotic 
medium property the need for which is more fundamental than anisotropy, is superluminal 



phase velocity^. In fact, all optical cloaking devices proposed so far, including the isotropic 
conformal cloak of Leonhardt et alJ ^^ * ^^ l, rely on the latter property, which in terms of re- 
fractive index can be stated as n < 1. True invisibility devices — the ones that can reduce 
the total scattering cross-section of an object situated in free space — preserve not only 
the amplitude but also the phase of all rays traversing through the device, which inevitably 
requires that at least somewhere in the cloak, there is a region with n < 1. The presence of 
such a region enables a curved ray to have an optical path equal to that of a straight ray in 
free space. 

Since the need for metamaterials in EM cloaking applications is mandated by this exotic 
property, one may realize that eliminating anisotropy from the implementation does not 
necessarily lead to a significant reduction in the design complexity, or a significant increase 
in performance. In fact, metamaterials and photonic crystals can be more easily and effi- 
ciently fabricated to have a strong anisotropy than to be perfectly isotropic. Metamaterials 
in particular are usually made of thin planar elements which provide strong response for one 
propagation direction only; eliminating the need for omnidirectional response allows one 
to stack them more densely and achieve higher strength of the response (in one direction 
only) and lower loss tangents. Thus, at least for two-dimensional optical cloaking applica- 
tions, reducing the range of transformations to conformal — and, consequently, the range 
of materials to isotropic — does not offer any substantial practical advantages. 

The real prospects offered by isotropic-medium designs open up in three dimensions, 
where anisotropic media are typically birefringent, that is, they support two normal 
electromagnetic modes with different dispersion relations. The coupling between the co- 
propagating modes is so difficult to manage that, in fact, it has been impeding the progress 
in three-dimensional transformation optics in the past five years. In the domain of elastody- 
namics (acoustics of elastic media), anisotropy leads to strong coupling^^ between longitu- 
dinal pressure (p) and transverse shear (s) waves, leading to issues similar to birefringence 
in electrodynamics. 

Unfortunately, conformal transformations in three dimensions are virtually non-existent, 
with one exception of the so-called sphere inversion transformation given by the famous 
Liouville's theorenP^. For that reason, the main focus of the recent conformal TO work 
has been on finding transformations that are approximately conformaP^'^ and thus lead 
to materials with anisotropy so weak that it can be neglected in the design, and then 



partially compensated by numerical optimizatioEp^. This branch of TO is now known as 
quasi-conformal TO (QCT0)IHEZ1. 

Eliminating anisotropy from the TWD designs, however, imposes serious — and quite 
fundamental — limitations to the performance of the resulting TWD devices. In particular, 
it is widely believed that a true, omnidirectional invisibility device cannot be obtained using 
only isotropic media, a belief articulated in the works of LeonhardiP^ and Greenleaf et al.'^. 
In this paper, we show that although cloaking devices based on conformal transformations 
are unlikely to become omnidirectional, for one specific angle of incidence they can indeed 
restore the amplitude and phase of the incident waves. 

This paper is organized as follows: first, we analyze the performance of the two- 
dimensional conformal cloak from Refs.'^^'^ using ray-tracing simulations that include the 
phase information, as well as full-wave simulations, and show that it can be made to per- 
form as a unidirectional invisibility device. Subsequently, we generalize that design to 
three-dimensions and propose an all-isotropic unidirectional cloak of a spherical shape for 
acoustic and electromagnetic applications. The EM (vector-wave) version of the proposed 
cloak is polarization-insensitive. 

II. ANALYSIS OF THE TWO-DIMENSIONAL CONFORMAL CLOAK 

Cloaking devices that arose from the TWD methodology are derived from coordinate 
transformations that compress certain regions of physical space into a zero-volume region 
in the virtual space. The two-dimensional variant of the cloak proposed by Pendry et al^^ 
shrinks a cylindrical volume isotropically to an infinitesimal diameter cylinder, a transfor- 
mation that requires anisotropic material properties but leads to a perfectly cylindrically 
symmetric distribution thereof. The inherently two-dimensional conformal cloak of Leon- 
hardt et al.'^SEl flattens a cylinder of radius a to a fiat sheet, which can be accomplished 
using a conformal map of an entire two-dimensional plane. Such conformal maps, known as 
Zhukovsky transform in hydrodynamics, have been applied to electrodynamic wave propa- 
gation problems long befor^^ the discovery of electromagnetic invisibility. 

Specifically, the transformation used here and in Refs.'^^'^ maps the exterior region r > a 
of a cylinder onto the entire plane with excluded line segment x G [—a, a]. The orientation of 
this segment, which we refer to as the principal axis of the conformal cloak (denoted x here), 
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FIG. 1. (color online). Two-dimensional conformal cloak proposed in Refs.'^^'^ (a) refractive 
index profile, with the boundary between superluminal (n < 1) and sub- luminal quadrants shown 
by green lines; the color shade inside the black circle (r = a) is cut off at refractive index values 
n = and n = 2, respectively; (b) refractive index as a function of the polar angle on the two 
sides of the discontinuity at r = a; (c) ray-tracing simulation showing trajectories and phase (color 
online) of the rays in the eikonal limit, assuming a plane wave incident along the x direction; (d) 
same as (c) with a plane wave propagating in the y direction. 



sets the direction of incidence for which this device may reduce visibihty of the cyhndrical 
cavity bounded by the circle r = a. In complex variables, this transformation can be written 

w{z) = z + a^ / z^ (1) 

where z = x + iy is the complex- valued representation of two-dimensional Cartesian coordi- 
nates, and w is the complex-valued coordinate in the virtual space. The practical result is 
that in the physical space, the refractive index distribution has to be 

^2 



ne{x,y) 



a^ — 2a^r^ cos 26 , , 

1 + -. -^ (2) 



{x + iy) 

where in the second expression we have used the polar coordinates (r, 0) oriented such that 
the positive x axis corresponds to = 0. This transformation is to be used only in the 
region exterior to the cylinder {r > a); what happens inside of it will be discussed below. 

The key feature of the refractive index distribution (pi) is that < rig < 2 for all r > a; 
the regions with superluminal phase velocity (n < 1) occupy two quadrants with roughly 
50% of the cloak volume, as shown by the isocontour n = 1 in Fig. flla). As discussed above, 
the presence of such regions is fundamental to the operation of any true invisibility device, 
as those regions enable optical path compensation. 

The main idea behind the device ^ is that a cylinder appears as a flat sheet, which — 
at least in the limit of geometrical optics — is invisible in the direction of the x-axis. Its low 
visibility for this particular angle of incidence is confirmed with the ray-tracing (Fig. fltc)) as 
well as full-wave (Fig. [2|a)) simulations performed with a commercial EM solver COMSOL. 
Of course, for any incidence angle other than = 0, the flattened sheet is not expected to be 
invisible, which makes the device given by the profile (pi) a unidirectional invisibility device. 

To increase the range of incidence angles for which the conformal mapping device (pi) can 
be invisible, it was proposed to extend the transformation onto the second Riemann sheet 
of the transformation (|T]) and ensure — by theoretical mean^^^JMj Q^j^gj; than transformation 
optics — that all rays emerging on the second Riemann sheet follow closed orbits. In 
order to implement this behavior, the interior of the circle r = a would have to be filled 
with inhomogeneous spatial refractive index distribution; for example, the Hooke's profile 
proposed in Refs.'^^'^ gives the following refractive index for r < a: 



n{x,y)=ne{x,y)\\\ —- , (3) 



where Ue is given by the same function ([2]), and w{z) — by the same equation Q. The 
complete index profile is shown in Fig. [lla). 

Notably, this index distribution is discontinuous across the boundary r = a, as illustrated 
by Fig. [l]^b). That such a discontinuity is unavoidable in the design methodology proposed 
in Refs .'23121 can be easily understood using the theory of functions of complex variables. One 
form of the uniqueness theorem for analytical functions of a complex variable states that 
an analytical function is uniquely determined by its values on a finite-length interval; the 
values of the function can be determined unambiguously using the process known as analytic 
continuation. Thus, if the refractive index n'{w) in the virtual space has to be continuous 
across the branch cut Hew = [—2a, 2a], it has to be given by the same analytic function 
on both Riemann sheets that are connected by that branch cut. Choosing n'{w) = 1 on 
one Riemann sheet forces one to either also pick n'{w) = 1 on the other sheet, or to have a 
discontinuity of refractive index. 

The behavior of the structure with refractive index dsl) filling the interior of the circle 
r = a in physical space is modeled using ray-tracing and full-wave simulations. To make a 
fair assessment of this structure, which is expected to work as a cloak only in the eikonal 
(short-wavelength) limit, we have chosen the free-space wavelength Aq = a/5 ^ a for the 
full-wave simulations, such that the conditions of geometrical optics are satisfied. Note that 
in the earlier assessment of this structure, the wavelength was chosen to exceed not only the 
diameter d of a cloaked object {d < a), but even the diameter of the entire cavity (2a); in 
that regime, the cavity and anything in it is poorly visible to begin with and thus its poor 
visibility is not sufficient evidence of the cloaking performance. 

While the ray-tracing simulations suggest that the full index profile (2|3 ) does indeed allow 



transmitted rays to return to their original trajectories after passing through the cloak, it is 
obvious that the transmitted rays leave the cavity (r < a) with incorrect field amplitude and 
phase^il. The phase delay is introduced to every ray that spends any amount of time inside 
the circle r = a, whereas the amplitude deviations are caused by the partial reflections on 
the refractive index discontinuity. These partial reflections are insignificant for the incidence 
angle = 0, at which the circle r = a is perfectly cloaked, but they become the limiting 
factor for all other incidence angles, at which the rays do hit the surface r = a from which 
they can get reflected. Note that the ray-tracing simulations presented here as well as in 
the original works^^EH completely disregard reflected rays, making an erroneous impression 
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FIG. 2. (color online). Two-dimensional confornial cloak proposed by Leonhardt et al! ^^ ^ I Full- 
wave simulations of the structure shown in Fig. fl} Angle of incidence: (a) 0, (b) vr/S, (c) vr/4, 
(d) 7r/2. TE polarization (E field out of plane) is assumed. Refractive index distribution is 
implemented using e^ = n^, with n given by Equations (2l3) and in-plane permeability /x = 1. 
Free-space wavelength Aq = a/5. 



that the transmission through this structure can be close to perfect. Strong reflectivity is 
revealed by the full- wave models shown in Fig. Il](c,d). 

To summarize, for any angle of incidence other than = 0, the device composed of index 
proflles ([2]) and ^ has a wide shadow, whose width at the 90° incidence is roughly twice 
the diameter of the circle (= 2a) and about four times the diameter of the cloaked cavity — 
the region of space not probed by any of the incident rays. This behavior is also observed 



at other angles of incidence, as evident from the full- wave simulations in Fig. [2] 

III. THREE-DIMENSIONAL ISOTROPIC-MEDIUM UNIDIRECTIONAL 
ACOUSTIC CLOAK 

The analysis presented in Section [IT] suggests that although isotropic metamaterial cloaks 
proposed in the earlier works can hardly be used as omnidirectional invisibility devices, for 
one angle of incidence they do have the potential to operate as cloaks. In two dimensions, 
they require the same type of hard-to-achieve response as their omnidirectional, anisotropic 
counterparts, namely refractive index in the — 1 range with very small imaginary part, 
and therefore are unlikely to beat the performance of the latter. However, extending these 
solutions into three dimensions would offer advantages such as elimination of birefringence, 
and, for acoustic applications, also reducing the coupling between longitudinal (pressure) 
and transverse (shear) waves. 

Three-dimensional extensions of two-dimensional TO and TA designs are difficult to 
come up with if the device has to be a precise implementation of a TO/TA prescription. 
However, in the so-called eikonal limit one only manages the trajectories and phases of 
the rays, allowing the impedance to vary almost arbitrarily in the device. Avoiding sharp 
discontinuities of the wave impedance allows one to achieve performance in eikonal-limit TO 
devices, including cloaks, that approaches that of the exact TO design. 

In the eikonal limit, a three-dimensional generalization of a TWD design can be obtained 
by a simple revolution of the refractive index distribution along some axis. For example, 
it is well-known that the refractive index profile of the omnidirectional 3D cloak of Pendry 
et al. is the same on any cross-section passing through its center as is the index profile 
of the 2D (cylindrical) cloak. Likewise, the ground-plane cloak that hides objects in the 
refiection mode can be quite easily generalized into 3D by rotation around a line normal 
to the ground plane. Here, we propose revolving the conformal cloak studied in Section [IT] 
around its principal axis, i.e. around the direction from which it can be actually invisible. 

As a result of the revolution, the cloaked cavity becomes spherical in shape, as seen from 
Fig. |3| Acoustic frequency- domain simulations based on the scalar Helmholtz equations can 
take advantage of the cylindrical symmetry of the problem using a standard axisymmetric 
acoustic solver^. The boundary condition on the surface of the spherical cavity is assumed 





(a) (b) 

FIG. 3. (color online). Three-dimensional unidirectional acoustic cloak consisting of isotropic 
medium: (a) pressure distribution on the cross-section of the cloak. Black lines show the streamlines 
of acoustic flux; acoustic wavelength at r ^ a is Aq = a/2; (b) three-dimensional picture of the cloak 
with a finite exterior radius Rcut = 4a; wavelength in ambient medium (at r > Rcut) Aq = a/2; the 
index transition is smoothed on a spatial scale AR = a/5. Acoustic wave impedance is assumed 
constant in both simulations; the speed of sound varies according to Cg = co/n, where n is prescribed 
by the revolution of profile pi) . 

to be Neumann {dp /On = 0), which corresponds to an acoustically hard scatterer, such 
as a solid object with a high bulk modulus (K) situated in a fluid much smaller bulk 
modulus {Kq <^ K). Very similar results (not shown) were obtained with Dirichlet boundary 
condition {p = 0), which is often used in acoustics to approximate an interface between 
two fluids. These two boundary conditions are also known as high- and low-impedance 
boundaries, respectively. As discussed in Ref.!^, for three-dimensional cloaks there is usually 
little difference between these two boundary conditions, since the surface is compressed and 
appears either as a point (for omnidirectional cloaks) or as a line (for the directional cloak 
described here) in the virtual space. In either case, point and line objects do not have a 
boundary condition associated with them, which explains why the structure is insensitive 
to the choice of the boundary condition in the physical space. 

Note however that allowing the waves to leak into the cavity by using a continuity bound- 
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ary condition still distorts the wave pattern and leads to additional scattering, which is 
evident from Fig. [2|a), where the interior of the cavity is filled with a refractive index and 
the surface r = a is modeled with a continuity boundary condition. As it turns out, in both 
2D and 3D the best performance is achieved when the surface of the cavity (r = a) is coated 
with a perfect reflector, as opposed to filling the cavity with the refractive index prescribed 
by©. 

So far we have been neglecting the fact that the refractive index distribution (pi) is filling 
the entire space, making this design impractical as an invisibility device. The simulation in 
Fig. |3](a), which shows nearly ideal cloaking, assumes that refractive index extends to infinite 
radius, and an artificial boundary (spherical PML) is introduced to mimic an infinitely 
extended cloak. A prior study^^ showed that simply cutting off the transformation at a 
finite radius Rcut results in a formation of a strong shadow generated by the refractive index 
mismatch 5n at r = Rcut- Since the index given by ^ scales roughly as 1 + a/r at r ^ a, 
the product 5nRcut does not converge to zero, and the effective scattering width of the cut- 
off cloak does not converge to zero in the limit Rcut ^ «• This type of scattering can be 
reduced, although not completely eliminated, by introducing a continuous transition from 
the index profile ([2]) to the homogeneous index of the ambient medium, which is normalized 
to unity in our calculations. The resulting simulation of a cloak with finite exterior radius 
is shown in Fig. [stb). 

To summarize, it is possible to achieve directional acoustic invisibility for a situation 
where the positions of the source and the detector relative to the cloaked object are known, 
using continuous distributions of isotropic density and isotropic bulk modulus. Such acoustic 
media are readily available as natural materials, random composites or mixtures, as well as 
ordered/structured metamaterials ^ * ^^ * ^^ *. Directional acoustic cloaks still require media 
with the speed of sound (cs = co/n, where cq is the speed of sound in ambient medium 
and n is the refractive index in the cloak) as high as possible^, and their performance is 
still limited by the highest achievable value of c^. We emphasize though that for limited- 
bandwidth cloaking applications only the phase velocity of sound needs to be high enough, 
and no statement about the group velocity is needed for such applications. In acoustic 
metamaterials, at a select frequency the phase velocity of sound can be engineered to be 
extremely high — much higher than the speed of sound estimated from the quasistatic bulk 
modulus and density of the materials involved, which opens up the possibility for nearly 
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FIG. 4. (color online). Three-dimensional polarization-insensitive electromagnetic cloak made of 
isotropic non-magnetic medium: (a) three-dimensional full-wave simulation with free-space wave- 
length Ao = 0.8a showing electric field component transverse to the propagation direction; simu- 
lation domain (excluding the PML) radius R = Aa] (b) 2.5-dimensional model with Aq = o/2 and 
cloak exterior radius Rcut = 6a; transition width Ai? = 0.3a. 

perfect acoustic cloaking^i^. 

IV. THREE-DIMENSIONAL POLARIZATION-INSENSITIVE 
ELECTROMAGNETIC CLOAK 



In the remainder of this paper, we apply the cylindrically-symmetric refractive index 



profile described in Section |III| to the problem of electromagnetic cloaking. As mentioned 
before, eliminating anisotropy from three-dimensional designs does away with birefringence, 
and allows a simple transition of TO solutions to their eikonal- limit. 

In order to stay within the validity of eikonal approximation, the cloaked object and 
the cloak itself must be many wavelengths in diameter. Modeling optically large structures 
with inhomogeneous refractive index tends to be a very computationally intensive task, and 
full-wave simulations of three-dimensional cloaks and other TO structures spanning many 
wavelengths are still not abundant in the literature. For a cylindrically-symmetric structure, 
one can reduce the computational domain to one quarter of the full domain by using two 
orthogonal symmetry planes intersecting along the axis of revolution. However, even with 
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that reduction, Finite Element simulations still often require prohibitively large resources 
and CPU time; the largest domain diameter to wavelength ratio we could efficiently model 
in this fashion was about 8, as shown in Fig. |4](a). 

Fig. llFa) reports a full-wave simulation of the directional cloak made of isotropic graded- 
index medium with the refractive index distribution obtained by revolution of ^ around 
the X axis. Our particular implementation of this refractive index in that simulation assumes 
a dielectric-only medium; similar results can be obtained with magnetic-only medium or a 
dielectric-and-magnetic medium whose refractive index satisfies the same equation. 

From the TO perspective, the described structure operates by compressing a spherical 
cavity covered with a perfect reflector (either PEC or PMC) to a zero-diameter, finite- 
length "needle". The needle is almost invisible from its tip in the geometrical optics limit, 
and it scatters little even when the wavelength is finite. Importantly, since the structure 
is rotationally-invariant with respect to the needle direction, the scattering pattern does 
not depend on the polarization of transverse EM waves incident upon it. As a result, this 
structure works as a directional cloak for an arbitrary polarization of incident light, including 
elliptic and partial polarizations. 

In order to characterize its scattering pattern deeply in the geometrical optics (eikonal) 
regime, we have developed a quasi-two-dimensional modeling method for vector EM fields, 
which can be referred to as 2. 5D modeling. The method is an extension of the well-known 
axisymmetric modeling technique used above for the acoustic wave cloak. The standard ax- 
isymmetric modeling assumes that all fields are independent of the azimuthal angle; there- 
fore, it cannot be applied to vector EM fields, whose transverse polarization breaks down 
the cylindrical symmetry even when the wave vector is aligned with the revolution axis 
of the structure. It is well known that EM scattering on a sphere (Mie scattering) has a 
distinct dipolar azimuthal dependence of the scattered field. The mathematical origin of 
this dipolar dependence can be seen by projecting the electric field vector of a plane EM 
wave onto the basis vectors of cylindrical coordinates: although the Cartesian component 
Er^ is independent of the azimuthal angle 0, the cylindrical components E^ = E^ cos (p and 
Er = E^ sin both have dipolar angular dependence. 

The 2.5D method differs from the usual axisymmetric modeling in that it includes a pre- 
scribed variation of the fields as a function of the azimuthal angle. In linear electrodynamics, 
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the fields can be decomposed into the Fourier series, 

E,,,,^(r,^,0) = 5^E(J,(r,^)e^-<^ (4) 

m 

and each cyhndrical harmonic propagates independently in the structure whose material 
properties are 0- independent. Note that in general all three field components are non- 
zero, and the full vector Helmholtz equation including the three unknown field components 
Ej^Ar^z) must be solved on the two-dimensional cross-section of the simulation domain 
for each m. After eliminating the 0-dependence, the equation depends on the azimuthal 
number m, and has to be solved for each cylindrical harmonic present in the incident field. 
Fortunately, if we are interested in modeling wave propagation at one angle of incidence 
(along the axis of revolution), only m = ±1 harmonics are non-zero in the incident field, so 
the equation needs to be solved only twice. In fact, the m = +1 and m = —1 harmonics are 
related by complex conjugation, which further reduces the computational load by a factor 
of two. 

The results for the directional three-dimensional cloak with exterior radius Rcut = 6a are 
shown in Fig. |4Fb), where the free-space wavelength is chosen to be Aq = a/2. The electric 
field pattern is very similar to the pressure pattern observed in the acoustic cloak shown 
in Fig. [stb). In this simulation, we have used incident electric field polarized in the = 
direction; due to rotational invariance of the refractive index distribution, it is sufficient to 
demonstrate cloaking for one polarization only, and the result extrapolates trivially to other 
linear polarizations, as well as their superpositions (elliptic polarizations). 

V. CONCLUSIONS 

In conclusion, we have analyzed the difficulties associated with the proposals to use 
isotropic media and isotropic metamaterials for cloaking applications. We have explained 
that achieving omnidirectional cloaking using isotropic-only media is extremely difficult even 
in two dimensions, where the design process can be guided by conformal transformation op- 
tics. Although no clear path to any such design has been identified so far, lifting the require- 
ment of all-angle invisibility allows one to arrive at useful two- and three-dimensional cloak- 
ing structures for acoustic and electromagnetic applications. Directional three-dimensional 
isotropic-medium cloaking designs eliminate the need to manage wave birefringence and 
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should therefore be substantially easier to make than their anisotropic counterparts. 
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